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We study the propagation of axial gravitational waves in Friedman universes. The evolution equa- 
tion is obtained in the Regge- Wheeler gauge. The gravitational waves obey the Huygens principle 
in the radiation dominated era, but in the matter dominated universe their propagation depends 
on their wavelengths, with the scale fixed essentially by the Hubble radius. Short waves practically 
satisfy the Huygens principle while long waves can backscatter off the curvature of a spacetime. 



I. INTRODUCTION. 
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Information about the topological type of the Universe can be obtained through the analysis of the Hubble relation 
[1]. One of the key and well known facts underlying derivation of the Hubble law is that electromagnetic waves 
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propagate along null cones, that is that they satisfy the Huygens principle in Friedman universes [2J. It is known that 
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gravitational waves do not obey the Huygens principle - they can backscatter on the curvature of the spacetime and a 
part of the radiation would come with a delay. The observation of the effects of the backscatter opens, in principle at 
least, a new local way to see the global topology in a way analogous to effects known in the Schwarzschild spacetime 
U3]-[5])_. 

The main goal of this paper is to explain whether and when these effects can manifest. We consider gravitational 
waves propagating in a fixed Friedman universe. The gravitational waves correspond to the so-called axial modes 
and they fulfill the Regge- Wheeler gauge condition [6] . The propagation equations become particularly simple in this 
gauge - they reduce to a single equation. This analysis with the Regge- Wheeler condition is novel in the cosmological 
literature, up to our knowledge. The form of the equation immediately implies that axial modes do satisfy the Huygens 
principle in the radiation epoch. On the other hand, the principle might well be broken in the matter-dominated 
epoch, depending on the wavelength of the gravitational radiation. The wavelength scale A is given by the average 
energy density p, A oc lA/p; under suitable conditions A is of the order of the Hubble radius - that is, of the radius 
■ of the visible part of the universe. The effects of backscattering are negligible in the short-range spectrum with 
wavelengths much smaller than A and they can possibly manifest for longer waves. That would mean that eventual 
effects cannot be observed by detectors of gravitational waves but they might be responsible for some superhorizon 
^ ' perturbations. 
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II. MAIN EQUATIONS. 



Let ds 2 = a 2 {—dry 2 + dr 2 + <r 2 d£l 2 ) be the line element of a Friedman cosmological model. Here r\ is the conformal 
time and the coordinate radius a = r, sinr and sinhr for the flat (k = 0), closed (k = 1 and open (fc = —1) universes, 
respectively. The axial perturbations of this metric in the Regge- Wheeler gauge [6] are defined as follows 

500 = h (ri, r) sin 9dgY w 

5i0 = hi(r],r) sin 9d e Yi ; (1) 

strictly saying, this a (^-independent subset of a larger class, but the general axial modes conform to the same 
dynamical equations as the ones displayed below. The linearized Einstein equations (for I > 2, applying arguments 
of [6] and [7]) are 

d n h Q = d r hi, (2) 
dlh x - d v d r h + 28^0— - 2d n (h^) + [I (I + 1) - 2] % = 0, (3) 
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d 2 h - d^drhx - 23^!— + 2^d r (a 2 ^) -1(1 + 1)!^+ 4kh = 0. (4) 
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This is a system of three equations on the two unknown functions ho and hi, but it is not over-determined, as we 
demonstrate in what follows. Acting on both sides of (4) by the operator <7 2 d v and using (2) in order to eliminate 
d v h in the obtained expression, one arrives eventually at 

d r F = Ak^dM - 2 (dlh x ) ad r a. (5) 

Here F is given 



F = a 1 (p 2 v hi - 8 2 r hi) +1(1 + 1) hi- 2o 2 % (hi^j . 
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Analogously, replacing d v ho by d r hi - thanks to (2) - in Eq. (3), one can obtain 

F = -2ad r ad r hi + 2hi. (7) 

Comparison of (5) and (7) allows one to conclude, that the two equations are equivalent - modulo terms in hi which 
are r-independent - if the compatibility condition 

-2d r hi [d r (ad r a) - 1 + 2ka 2 ] = 0. (8) 

is satisfied. It is a straighforward exercise to check out that the expression within the square bracket vanishes for all 
Friedman models. If one assumes in addition compact support of initial data, then the above mentioned ambiguity 
disappears and one can conclude that the two equations, (5) and (7), are equivalent. The whole dynamical content 
of the evolution is described by Eqs. (6) and (7), 

a 2 (d 2 hi - d 2 hi) +1(1 + 1) hi- 2<r 2 d ri (hi^^j = -2od r od r hi + 2h x . (9) 

Insert now Q = hi/(aa) into equation (9); this yields 

d 2 Q-d 2 r Q+ l -^^-Q-kQ-^Q = 0. (10) 



One can show, using Friedman equations, that in the radiation epoch d 2 a = — ka while in the dust epoch d 2 a = 
— ka + 47r/?a 3 /3. Thence the gravitational wave evolves during radiation era according to the equation 

d 2 Q-d 2 Q+ l -V±±±Q = 0. (11) 
During the matter (and, in particular, dust) dominated epoch, in turn, the gravitational waves satisfy 

d 2 v Q-d 2 Q+ l ^Q-^Q = 0. (12) 



III. ANALYSIS. 



During the radiation epoch the Huygens principle holds true. Indeed, the form of the equation (11) is exactly that 
of the electromagnetic fields in Friedman spacetimes. Therefore its general solution has the form (see, for instance, 
[8]) 

Mr,V) =° l d r -d r ----d r (t±^) (13) 

a a a 

V v ' 

/ 

where the free functions / and g depend on the combinations r — r\ or r + r\, respectively. Clearly, from this form 
of the solution, compact pulses of radiation have to move within space-time regions bounded by a pair of null cones 
i] = +r + const. In particular shocks have to move along null cones. 

Now we shall discuss the evolution of gravitational waves in the matter-dominated era. Assume initial data of 
compact support a(r) < cr(r ) at the beginning 770 of this epoch. Thus at a time rj the radius of the support of Q is 
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bounded from above by a(ro+r]); assume in addition a(r + r/) < tt/2 and k/a << H in the case of the k — 1 universe. 
Let R(r)) = a(rj)(j(r + 77) be the related areal radius. Define the critical-length parameter X(rj) = y/3l(l + l)/(4np). 
Let us remark that, in the case of the flat universe, A (77) is of the order of the quantity y/2l (I + 1)/H (where H is the 
Hubble constant). This is related to the Hubble radius and for low valus of I is of the order of the radius of the visible 
part of the universe. In the case of closed and open universes X(rj) can be bounded from above or below, respectively, 
by y/2l (I + Tj/H - again, for small I, being close to the Hubble radius 1/H. One can write, taking into account that 
the product p(rj)a 3 (rj) is constant during the matter dominated era, a sequence of inequalities relating the ratio N(rf) 
of the two radii, taken at succesive moments of conformal time, 

« 2 (*)«fe) =J v M(1 + 2(1 , 2)) . (14) 



A 2 (r/ 2 ) a(»?i) 

2 3 i o + Vi d /lAOTYinlArn n ■ 

a(ri 



Here z(l, 2) = is the cosmological redshift and 770 < 771 < rj 2 - 



The ratio of the last two terms of Eq. (12) is bounded from above by 4k pa 2 (r\)o 2 (r + 77)/ (31(1 + 1)). This can 
be expressed as N(r]) = R 2 (77)/ A 2 (77). Clearly, if AT (772) (1 + z(l, 2)) << 1, then N(rji) « 1. That suggests that 
for 77 G (770,772) the evolution is well described by the equation (11) and the solution is close to that given by (13). 
Therefore there is no backscatter. And conversely, if the above product is at least of the order of unity, then there can 
exist a backscattered part, with characteristics possibly allowing the identification of the type of the universe. This 
heuristic reasoning can be made rigorous, as outlined below. 

One can easily check that the energy 



E{Q) ^IJ dr 



(d v Q) 2 + (d r Q?+ l ^±±±Q 2 



(15) 



is conserved by the evolution equation (12). Split the sought solution of (12) as Q — Qo + x, where Qo satisfies Eq. 
(11) and initially Q = Qo, d v Q — <9 r) Qo- Thence x solves the equation 

dx-d r x-\ 5 — x= — - — (x + Qo) (16) 

a o 

and satisfies homogeneous initial conditions x = d v x = 0. The energy E(x) of the field x satisfies the equation 

±E(x) = J™ dr^-(Q + x)d ri x. (17) 

One notices that a (ro + 77) > a (r) for all points r lying within the support of Qo- 
Therefore equation (17) implies the inequality 



i E(x) s ijllr^f *^ IWo+ * )fW (18) 

Now, making use of the Schwarz inequality and bounding the relevant factors on the the right hand side of (18) by 
energies E(x) and E(Qq), dividing both sides by \J E(x) and rearranging the resulting inequality, one obtains 

e a(„) JL (e-WyiM) < ^VEW), (19) 

where 

In the case of flat and open universes one has An pa 2 / 3 < H 2 a 2 /2, while for the closed universe there is an approximate 
equality 47rpa 2 /3 w H 2 a 2 /2, under the condition imposed earlier. Notice now that dr/Ha 2 — da; change the integration 
variable and observe that a (r + 77) and H (rj) are increasing with 77, in order to get the estimate 
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a < a (r + v) (« (v) ~ a fa,)) - |M= . (21) 

2V t (t + 1) 



This can be written concisely in terms of the areal radius and the critical parameter, 



Finally, one arrives at the inequality 



a<^ . (22) 



S(x,r?) < J B(g )(e /3 -l) 2 . (23) 



Similar estimates can be obtained for the energy E (d r x) and of higher derivatives of x. Employing these, one can 
use standard Sobolev estimates [9] in order to show that the field x (and its derivatives, if needed) can be pointwisc 
small, under the condition that initial data at the time r] = have compact support and that R (rj) / (A (rj)) « 1. 

Let us remark that we omit boundary terms coming from r = 0, without discussing subtleties related with the 
behaviour of solutions at r = 0, when integrating by parts. This is justified by the fact that "good" initial data of 
our 1 + 1 evolution problem - i. e., those inherited from the initial 1 + 3 formulation of the evolution of gravitational 
waves - take care of themselves at r = and they give no contribution in relevant boundary terms that might appear 
above. 

The compactness condition R(r,rj) « A (rj) means that the areal size of the radiation pulse is small comparing 
to the critical radius A (rj). This in turn implies, using the Fourier transforms [10], that in this (no-backscatter) case 
the spectrum is dominated by high frequencies. The compelling explanation of the observed correlation (dominant 
high frequency)-(no backscatter) is that short length-waves do not undergo backscatter. If the above condition is 
not satisfied, then we do not have analytic arguments that would necessitate strong backscatter, but a preliminary 
numerical analysis suggests that this indeed is the case. 

In conclusion, short gravitational waves - with the adjective "short" meaning a wavelength much smaller than the 
Hubble radius for the k = 0, — 1 universes - propagate without backscatter. Longer gravitational waves, in particular 
those with the wavelengths being of the order of the radius of the horizon (for k = 0, — 1 Friedman models, and also 
for the k = 1 model, under additional conditions) are expected to backscatter during the matter dominated era. This 
effect can manifest most strongly in the superhorizon scale. 
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